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An experimented study of plasmachemical reaction involving CH 4 and N 2 molecules 
in rf discharge was studied in order to know the effect of vibrational excitation of 
N 2 molecules. When the relative nitrogen concentration was greater than 0.8, the 
main product of CH 4 decomposition was HCN, and the rate of methane decompo- 
sition at this condition was faster than that one in pure methane. These results 
could be confirmed through the mass spectroscopic method. The reason for these 
results is the vibrational energy of N 2 excited by rf discharge. The chain reaction 
mechanisms of producing HCN by vibrational excitation of N 2 were examined 
closely through numerical simulation. The rate-controlling step was the dissociation 
reaction of excited nitrogen molecule to the atomic nitrogen, so the process of HCN 
synthesis was limited by the value of reaction constant, k N . 
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1. INTRODUCTION 

Studies of chemical reactions in non-equilibrium molecular plasma at 
elevated pressures have been closely related to the progress of plasmachem- 
istry, hydrogen power engineering, waste-handling of natural gases, cleaning 
of an environment, etc. The energy efficiency of non-equilibrium plasma- 
chemical process depends on the set of channels it flows, i.e., on the mechan- 
ism of the process. It has been known that the vibrational excitation of 
molecules essentially accelerates endothermic chemical reactions. 11 ’ How- 
ever, it is not always possible to excite the required vibrational mode of 
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molecules selectively by an electric discharge. In our previous work/ 2 ’ 3 ’ we 
investigated the decomposition of pure methane and carbon dioxide in a 
radio-frequency discharge. It was shown that the dissociation of these mol- 
ecules was due to the excitation of electronic states. The plasmachemical 
reactions in nitrogen mixtures were examined in order to analyze the effect 
of vibrational excitation on the reactions involving methane. 

Molecular nitrogen has a large effective cross section of the vibrational 
levels excited by electron impact (3 x 10 _l6 cm 2 ), a small effective cross sec- 
tion of vibrational relaxation (3 x 1 O' 24 cm 2 ) and a small factor of the 
vibrational energy loss on the surfaces. For glass, quartz, stainless steel and 
copper, this factor for the accommodation of vibrational energy loss is equal 
to about lOV" In other words, N 2 molecules are excited in discharge very 
easily and act as a reservoir of the vibrational energy. 

2. EXPERIMENTAL 

We investigated the plasmachemical reactions involving CH 4 and N 2 
molecules in radio-frequency discharge (v= 13.56 MHz) by a mass spectro- 
scopic method. These reactions took place in discharge in the gas mixtures 
of CH 4 and N 2 . We used a special type of capacitive discharge. A similar 
discharge system was applied at first for the design of C0 2 lasers by 
Yatscnko' 4 ' and was later used in this experiment for plasmachemical 
purposes/ 2 ’ 

The schematic drawing of the experimental setup was shown in Fig. 1. 
Plasmachemical reactor consists of a long pyrex (or quartz) tube. Four cop- 
per wires were located on the outside tube and were used as electrodes. The 
diameter of each wire was d/ 10, where d is the inner diameter of the reactor. 
Any two of these were connected with power supply and the other two were 
connected to earth. The reactor was made of Pyrex glass with an internal 
diameter of 12 mm and a total length of 700 mm, and the plasma zone had 
500 mm length. More detailed descriptions of the plasmachemical reactor 
and all experimental equipment were described in our previous work/ 2 ’ The 
main peculiarity of these reactors was the small sizes of the electrode 
sheathes. As a result almost all volume of the discharge tube was filled with 
positive column plasma/ 4 ’ The pressure of gas mixture was changed from 5 
to 60 torr. 

The radio-frequency generator with a matching network delivered an 
output power from 0 to 300 W. The magnitude of reflected power did not 
exceed 2% from the delivered one. The maximum of unique power for the 
reactor was about 7.2W/cm 3 . While measuring discharge input power, we 
ignored the energy loss through radiation and, furthermore, we suggested 
that all input power was absorbed by positive column plasma. 
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Fig. 1 . Schematic drawing of the experimental setup. 



CH 4 and N 2 with 99.9% purity were used. Quadrapole mass spec- 
trometer (Balzers, QMS 200) with Quadstar 421 software was used for 
qualitative and quantitative analysis of the gas mixtures. Mass spectrometer 
was connected to the post discharge zone. Gas mixture in this zone was 
maintained at room temperature. Before measurements, we carried out a 
calibration of the mass spectrometer with the data based on the mass spec- 
trum of the binary mixtures. CH 4 , C 2 H 6 , C 2 H 4 , C 3 H 8 , and Ar gases with 
99.9% purity were used for the calibration. 

Some expressions in our previous world 2 ’ were used to define the con- 
versions of initial reactants and molecular flows of reactants and products 
investigated. The residence time was considered for the change of the flow- 
rate by chemical reactions. 

3. RESULTS AND DISCUSSION 

Let us consider plasmachemical reactions in discharge in mixture 
of CH 4 and N 2 . The effect of relative nitrogen concentration 
(/3 N2 = [N 2 ]°/[N 2 ]° + [CH 4 ]°, here [N 2 ]° is the initial concentration of the 
nitrogen molecules and [CH 4 ]° is the initial one of the methane molecules) 
was investigated on plasmachemical processes. The mass spectra analysis 
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showed that the influence of nitrogen was minor for /3 N2 < 0.65. The situ- 
ation was very similar to the discharge in pure methane/ 2 ’ At low input 
power of 120 W ethane and hydrogen were the main products. As the input 
power was increased, the unsaturated groups of C 2 and C 3 began to form. 
At j 3 N2 >0.8 the situation was quite different. The main products of CH 4 
decomposition were HCN and H 2 . No other substances were detected in 
noticeable amounts. 

At /3 N2 >0.8 the main plasmachemical process in the discharge is as 
follows: 

CH 4 + ^N 2 ^HCN + |H 2 (1) 

As an illustration of this process, the dependencies of the methane con- 
version and the ratio of Fh 2 /^ch 4 at relative nitrogen concentration are pre- 
sented in Figs. 2 and 3 under various input powers. Here, Th 2 is the flow 
rate of the molecular hydrogen in the post-plasma zone, /vn 4 is that of the 
methane molecules in the chemical reaction zone. Initial conditions were that 
total pressure PflO) was 23 torr and total flow rate V° was 55 cm 3 / min. It 
should be noted that the methane conversion increased as the value of /3 N2 
increased. At /3 N2 = 0.9 and IF>200W almost all methane was converted 
into HCN and H 2 (Z CH4 > 0.9, where Z CH4 is the conversion of methane). 

The ratio of Fh 2 /Fch 4 , was almost constant at IF>200W, and the 
value was equal to 1.5 as shown in Fig. 3. At 120 W of input power and 




Fig. 2. The dependencies of methane conversion Z on relative nitrogen concentration for 3 
values of input power. The initial conditions were: discharge in mixture CH 4 -N 2 , total pressure 
is /MO) = 23 torr, and total flowrate F <0) = 55 cm 3 /min. 
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P=[n 2 ]/([n 2 ]+[chj) 

Fig. 3. The dependency of the ratio FuJ Fch 4 on relative nitrogen concentration for 3 values 
of input power. The initial conditions were: discharge in mixture CH 4 -N 2 , total pressure 
PdO) = 23 torr, and total flowrate V (0) = 55 cm 3 /min. 

/3 N2 <0.65 the value of Fu 2 /Fch 4 was equal to about 0.8. At /3 N2 >0.8 this 
value was equal to 1.5. In the region of 0.65 < /3 Nz < 0.8, it showed a drastic 
change of Fh 2 /Tch 4 from 0.8 to 1.5. The reaction mechanisms changed 
in this region. If /3 N2 <0.65 the nitrogen molecules were not involved in 
plasmachemical reactions, but at /3 N2 >0.8 the chain reaction occurred as 
follows:* 1 ’ 



CH 4 + N- 


A-hcn 

> 


hcn + h 2 + h, 


AH =-0.51 eV 


(2) 


H + Nf- 


A:N 

> 


• NH + N, 


£„ = AH = 6.54 eV 


(3) 


NH + NH- 


fc N 2 

) 


•n 2 + h 2 , 


AH = -7.8 eV 


(4) 



where AH is the standard enthalpy of the reaction, and E a is the activation 
energy of the reaction. It can be shown that the next combined rate equa- 
tions describe this mechanism (2)— (4) : 

_ r/[CH 4 ] = ^HCN [N][CH4 j ( 5 ) 

dt 

k K C t [N 2 ] - (/c n [N 2 ] + /i" cn [CH 4 ])[N] 
dt 



( 6 ) 
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= k N 2 [NH f - k N (C t - [N])[N 2 ] (7) 

dt 

= /r N (Ci - [N])[N 2 ] - 2 /t N2 [NH] 2 (8) 

dt 



where the values in the square bracket are the concentration of the relevant 
substances, k HCN , k N , and k Nl are the rate constants of reactions (2)-(4), C i 
is a constant value. (It can be derived from reactions (2)-(4) that C i = 
[N] + [H].) 

It is obvious, that 

[HCN] = [CH 4 ]° - [CH 4 ] (9) 

where [CH 4 ]° is the initial concentration of methane molecules (that is, 
[CH 4 ]° is the concentration in the predischarge zone), and [CH 4 ] is the cur- 
rent concentration. 

Let us consider the peculiarities of the processes (2)-(4). It is easy to 
understand that the synthesis is limited by endothermic reaction (3). This 
reaction is stimulated by vibrational excitation of nitrogen molecules quite 
well. It is a reasonable assumption that k N «k HC1<l ~k™ 2 . 

The effect of atomic-nitrogen concentration is our primary concern. 
From Eq. (6), we can find 

k N Ci[N 2 ] - C 2 exp(-(7c N [N 2 ] + /c ncN [CH 4 ])r) 



[N] = 



/r N [N 2 ] + At" cn [CH 4 ] 



( 10 ) 



Where C 2 is a constant value defined from the initial condition of 
[N (t = 0)] = [N]°. From Eq. (10), it may be seen that under the condition of 

t>h = (/c n [N 2 ] + /c ntN [CH 4 ]r' 



the atomic-nitrogen concentration will be 

^ N [N 2 ]Ci 



[N]i = 



jCH 4 ] + /t n [N 2 



at any initial conditions. If 

k l l(N [CH 4 ]° » /r N [N 2 ]°, ?, = (/r ncN [CH 4 ]r' and [N], = 



[N 2 



r CN [ch 4 ] 



(ii) 



( 12 ) 



Cl 



Let us consider now the rate equation (5) for the methane concen- 
tration. From Eq. (5), the time scale of the substantial methane concen- 
tration change is t~(Ci& hcn ) _1 , where t is the residence time of the reactor. 
Notice that it is estimation for the minimum time. 
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In the conditions of the molecular gas discharge plasma, the concen- 
trations of atomic species are much less than the initial concentrations of 
the molecular species. 1 (6) Hence, the residence time is greater than t t (t»C), 
the concentrations of [CH 4 ] and [N 2 ] are almost the same during the period 
of time from 0 to . It means that we can substitute the initial values of 
[CH 4 ]° and [N 2 ]° into Eq. (12) and the value of [N] t from Eq. (12) into Eq. 
(5), thus we have the next equation, 



4CH 4 ] [n 2 ]° 

— K A C^i dt 

[CH 4 ] [CH 4 ]° 



(13) 



From the Eq. (13), the methane concentration decreases by the following 
equation, 



[CH 4 ] = [CH 4 ] u exp -jrCi 



[N 2 



[CH 4 ] u 



(14) 



with the characteristic time, 



1 [CH 4 ]° 

/c N C, [N 2 ]° 



(15) 



It is necessary to point out that t 2 is the time for the HCN production 
(Eq. (2)). 

Thus, as mentioned in the above relations, the plasmachemical process 
of HCN synthesis was independent of the initial value of [N]°. During the 
time period of order of fi~(/r HCN [CH 4 ]°) _1 , the concentration of atomic 
nitrogen defined by Eq. (12) is used in the system under investigation. The 
densities of [CH 4 ]° and [N 2 ]° do not change practically during this time. 
Then during the time period of order 



( 1 [CH 4 ]°\ 

U N c, [n 2 ]° ! 



»t 1 



the concentration of methane molecules decreases noticeably. (In this time 
the noticeable amount of [HCN] is produced.) It is obvious that for r »t 2 , 
when the concentration of CH 4 decreases significantly, the value of [N] will 
be equal to Ci (see Eq. (12)) and the concentration of [CH 4 ] will decrease 
by the following relation, 



[CHJ-exp^n^CdO 

with the characteristic time of T- (/r HCN Ci)“ 1 . 



( 16 ) 
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Fig. 4. The time on stream of N 2 , CH 4 , NH, and N concentrations. The initial conditions 
were: [CH 4 f = 0.091, [N] 0 ' = 1(T 6 , [N 2 f = 0.091, and [NH] 0 ' = 1(T 7 . 



For the examination of the above estimations, we made a numerical 
modeling of the process described by Eqs. (5)— (8). To make a simple analy- 
sis more easy we introduce the new variables for the concentrations and the 
rate constants, []'=[ ]/[CH 4 ]° + [N 2 ]°. The unit of [ ]' is dimensionless. 

In order to solve these equations, we used the Runge-Kutta method. 
The initial densities were [CH 4 ] 0 ' = 0.091, [N 2 ] 0/ = 0.91, [NH] 0 ' = 0, and the 
concentration of atomic nitrogen [N] 0 ' was changed. All results are pre- 
sented in a graphical form with dependence of ln[ ]' on time. In the condition 
under investigation ( P R = 23 torr, V° = 55 cm 3 / min, T R = 800 K) the mean 
residence time for molecules in plasma was about t- 0.5 sec, which was 
defined as the characteristic time scale. 

Figures 4 and 5 show the dependencies of results on initial value of 
[N] 0 '. For Fig. 4 the value of the nitrogen atom concentration was [N] 0 ' = 
10 -6 , and for Fig. 5 that one [N] 0 ' = C i = 1 0~\ From these figures two results 
were obtained. 

Firstly, practically there is no influence of the initial density of the 
nitrogen atoms on the time dependency of [CH 4 ]' (or on the time depen- 
dency of [HCN]). Secondly, it is possible to use the simple estimation for 
the times t, = 1.4 x 10“ 3 sec and t 2 = 2.5 sec from Eqs. (1 1) and (15). 
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Fig. 5. The time on stream of N 2 , CH 4 , NH, and N concentrations to show the influence of [N] 0 '. 
The initial conditions were: [CH 4 ] 0 ' = 0.091, [N] 0 ' = Ci = 10 3 , [N 2 ] 0 ' = 0.91, and [NH] 0 ' = 1(T 7 . 



Figure 6 shows the effect of the rate constants of k HCN and /r Nz with 
the reaction time. In this figure the values of /c nCN and k^ 2 were greater 
than those in Figs. 4 and 5 by a factor of 10. These figures showed that the 
time of t 2 in Fig. 6 was approximately equal to those in Figs. 4 and 5. 
When the value of [N]' was equal to C i = 0.001, the rate of the methane 
decomposition increased more sharply in Fig. 6 than Figs. 4 and 5. But this 
fact had no practical importance, because the methane concentration had 
decreased at this moment by two orders of magnitude (i.e., the methane was 
almost decomposed) and most of HCN had been produced. Hence the pro- 
cess of HCN synthesis is not strongly affected by the values of k HCN and 
k N \ 

Figure 7 demonstrates the effect of the value of /r N as time goes by. In 
this figure k s was increased by a factor of 4 in comparison with previous 
figures. The time t 2 = 0.63 decreased accordingly by a factor of 4, and the 
production of HCN was accelerated noticeably. 

Up to this point, we considered that A' HCN = k Nz . Figures 8 and 9 dem- 
onstrate that this assumption is not crucial. In conditions under investi- 
gation, when the value of [N 2 ] is noticeably greater than the value of [CH 4 ], 
the value of k^ 2 has no practical effect on time evolution of the methane 
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Fig. 6. The time on stream of N 2 , CH 4 , NH, and N concentrations to show the influence of 
values of /c HCN and A: N2 . The initial conditions were: [CH 4 ] 0 ' = 0.091, [N] 0 ' = 10 3 , [N 2 ] 0 ' = 0.91 
and [NH] 0 ' = Hr 7 . 



decomposition. k Nl - 180,000 sec -1 in Fig. 8 and k Nl - 8000 sec -1 in Fig. 9. 
Nevertheless the time of changing [CH 4 ] and [N] are the same for these 
figures. The value of /c Nz has an effect only on the value of the intermediate 
product density, [NH], but has no noticeable effect on the rate of 
decomposition. 

Thus the results of the numerical modeling supported the validity of the 
estimations which were made on the basis of the simplified consideration. 

In our previous work, <2) we obtained the expression for describing 
methane decomposition process in discharge with the pure methane. In that 
case, the methane-concentration change was described as follows 

[CH 4 ] = [CH 4 ]° Qxp{-n e (v e a e di SS )t} (17) 

where t is the residence time. The frequency of collisions for methane mole- 
cules with electrons is n e (v e a e diss ) where n e is the density of electrons, v e is 
the speed of electrons and Od iss is an effective cross section for dissociation 
by direct electron impact. 

A comparison between Eq. (17) and Eq. (14) shows that the methane 
concentration decays exponentially in discharge with pure methane and with 
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Fig. 7. The time on stream of N 2 , CH 4 , NH, and N concentrations to show the influence of 
values of £ N The initial conditions were: [CH 4 ] 0 ' = 0.091, [N] 0 ' = Ci = 1(T 3 , [N 2 ] 0 ' = 0.91, and 
[NH] 0 ' = 10~ 7 . 



the mixture of CH 4 and N 2 . If 

[Cri 4 J 

the nitrogen molecules are not involved in plasmachemical process. If 

, , , /N ^ [N 2 ]° 

n e (v e o diss ) < k C, 

|CH 4 J 

the mechanisms of (2)— (3) are signihcant and nitrogen molecules are 
involved in chemical reactions. 

If /3 N2 < 0.65 the nitrogen molecules were not involved in plasma- 
chemical reactions. The production of nitrogen atoms was too small. But 
at /3 N2 >0.8 the chain reaction occurred. The reaction mechanisms 
changed in the region of 0.65 < /3 N2 < 0.8. The discussed mechanism was 
not influenced by the total concentration of nitrogen and hydrogen atoms 
(Ci = [N] + [H]). The presence of initial nitrogen and hydrogen atoms in 
plasma could be determined by the electron impact dissociation of N 2 
and H 2 . But it should be mentioned that the electron impact dissociation 
of N 2 in rf discharge could not play the main role in the production of 




170 



Savinov el al. 




Fig. 8. The time on stream of N 2 , CH 4 , NH, and N concentrations to show the influence of 
values of £ HCN and k N2 . The initial conditions were: [CH 4 ] 0 ' = 0.091, [N] 0 ' = C,= 10 ~ 3 , [N 2 ] 0 ' = 
0.91, and [NH] 0 ' = Hr 7 . 



N atoms at the condition under investigation. The efficiency of electron 
impact dissociation depended on the parameters of gas discharge. The 
reduced electric field strength E/N (E is the longitudinal electrical field 
strength and N is the density of neutral plasma gas particles) was one 
of major parameters. This parameter determined the mean energy of 
electrons in plasma and, accordingly, the energy balance of different 
degrees of freedom of the atomic and the molecular components in 
plasma. But under changing /J N2 from 0.65 to 0.8 E/N changed very 
lightly and the efficiency of electron impact dissociation could not rise 
noticeably in this region. 

It is necessary to point out that the dissociation energy of N 2 is higher 
than the bond dissociation energy of CH 3 -H by a factor of 2. Nevertheless 
under favorable conditions the N 2 is excited easily up to high vibrational 
levels. The presence of many vibrational excited molecules noticeably 
increases the value of /c N . When 

«,(v,aU<* N C,-^L 

[th 4 j 

the chain mechanism of (2)— (4) began to dominate and the synthesis of 
HCN took place in the plasma. Because & n Ci([N 2 ]°/[CH 4 ] 0 ) was noticeably 
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Fig. 9. The time on stream of N 2 , CH 4 , NH, and N concentrations to show the influence of 
values of /: HCN and k Nl . The initial conditions were: [CH 4 ] 0 ' = 0.091, [N] 0 ' = Ci = 10“ 3 , [N 2 ] 0 ' = 
0.91, and [NH]°'= 1(T 7 . 



larger than n e (v e o e d i ss ) in this condition, the transition to this mode of dis- 
charge operation had to be accompanied by the increase of the methane 
conversion. 



4. CONCLUSIONS 

An experimental study of plasmachemical reaction involving CH 4 mol- 
ecules with N 2 in rf discharge was investigated with a mass spectroscopic 
method. When the relative nitrogen concentration was greater than 0.8, the 
main product of CH 4 decomposition was HCN. However, the other con- 
ditions, especially the relative N 2 concentration was less than 0.6, the reac- 
tion mechanism was almost the same as the decomposition of pure CH 4 . 

The chain reaction mechanisms of producing HCN by vibrational exci- 
tation of N 2 were examined closely through experiments and numerical 
simulation. The rate-controlling step was the dissociation reaction of excited 
nitrogen molecule to the atomioc nitrogen. As a result, practically there was 
no influence of the initial concentration of the nitrogen atoms on the time 
dependency of CH 4 or HCN concentration. 
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When the value of rate constant, k N , increased, the reaction time 
decreased accordingly and the production of HCN was accelerated notice- 
ably. Hence the process of HCN synthesis was induced which was limited 
by the value of k K . The values of k HCN and k Nz affected only the concen- 
tration of NH, but have no noticeable effect on the rate of decomposition. 

In the frame of this work we had no possibility of defining the values 
of Ci and k N . By using Eq. (14) and the value of mean residence time for 
molecules in plasma (x- 0.5 sec) we could obtain information about the 
product of C|/i N ~0.5 s _1 . In the future it will be interesting to make special 
experiments to measure the vibrational temperature of N 2 and to estimate 
the atomic concentration in the plasma. It enables us to obtain new infor- 
mation about the process under investigation. 
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NOMENCLATURE 

Ci = constant, [H] + [N] 

E a = activation energy, eV 
F - flow of molecular spices, min -1 
AH = standard enthalpy for reaction, eV 
k = rate constant, sec' 1 
L = reactor length, cm 
n e = electron density, cm -3 
P = pressure, torr 

S = cross section of plasmachemical reactor, cm 2 

T = temperature, K 

V = volumetric flowrate, cm 3 min -1 

v e - electron velocity, cm sec -1 

W = specific input power, W cm -1 

Z, = conversion of i component 

Greek Letters 

Pi = relative concentration of i component 
8 r = factor for calibrating flowrate change by chemical reaction 
£ = specific energy, J 

CTdiss = effective cross section of collision for dissociation by direct elec- 
tron inpact, cm 2 
t = residence time, sec 
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Superscripts 

HCN = hydrogen cyanide 
N = atomic nitrogen 
N 2 = nitrogen molecule 
Pr = product 
R = reaction 
0 = initial state 

Subscripts 

e - electron 
R = reaction 
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